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Abstract - All present ferroelectric transistors "have been made on the micrometer scale. 

Existing models of these devices do. not take into account effects of nanoscale 
ferroelectric transistors. Understanding the cb^acferistics of these nanoscale devices is 
important in developing a strategy for building and using future devices. This paper 
takes an existing microscale ferroelectric field effect transistor (FFET) model and adds 
effects that become important at a nanosca|e,.leyel, including electron velocity saturation 
and direct tunneling. The new model analyzed. FFETs ranging in length from 40,000 
nanometers to 4 nanometers and ferroelectric, .thickness from 200 nanometers to 1 
nanometer. The results show that FFETs can operate on the nanoscale but have some 
undesirable characteristics at very small dimensions. 


Introduction 


The purpose of this paper is to develop a modelithat can accurately predict the 
characteristics of ferroelectric field effect transistors. (FFET) for a range of sizes of 
channel lengths from tens of microns to a fgW UaUpmeters. Models exist [1] that can 
accurately characterize FFETs in the micron' scale size. These can be verified to be 
accurate by comparison with data from actual, FFETs, Recent research [2] [3] has 
developed equations and models of very splqlljransistors (non-ferroelectric). This paper 
takes an existing FFET model and incorpqrafcs ^jfifuiets from having a very short channel 
length and very thin ferroelectric films over tbe channel. The model simulates only 
classical FET layout design of a transistor- and .ctaes not attempt to simulate single 
molecule or vortex layouts. The basis forsthd empirical data is a PZT FFET w with a 
channel length of approximately 40 microns. The model was run for channel lengths of 
40,000 nm, 4000 nm, 400 nm, 40 nm, and 4 nm-uTo reduce complexity all W/L ratios 
were 1. The ferroelectric film thicknesses ranged from 200 nm to 3 nm. For thicknesses 
of less than 3 nm, the model did not produce, ^q,f ^fesp|ts.' This may be due to instability 
in the model or actual physical problems dpeTqThe short distances and high electric 
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fields. Ferroelectric films of less than 2 nm would be difficult to produce and the 
perovoskite crystal structure has physical size limitation in this range. 


II. Basis ofthe Model 


A. Partitioned Ferroelectric Model 

The major difference between the Ferroelectric Field Effect Transistor and the normal 
MOSFET is the replacement of the oxide layer with a ferroelectric layer above the 
channel of the transistor. This layer develops charge much differently than the oxide 
layer due to polarization of the material. These charges affect the surface potential within 
the channel and change the threshold voltage of the transistor. Several researchers have 
studied these effects and have developed mpdels o|how they work. Wu 15J developed a 
mathematical model that accounts for the additional charges from the ferroelectric 
material at the semiconductor-ferroelectric interface. Since the threshold voltage for a 
MOSFET is a function of the total charge ' ^t the $emiconductor-insulator interface, the 
threshold voltage is: ' y 
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Where Ps is the ferroelectric polarization, Ap is the ferroelectric material area, and Cd is 
the insulator capacitance (or ferroelectric capacitance, Cf). This new Cd is calculated 
using the following equation. 
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Where, C ox is the capacitance of the oxidel^yer, (4 c ; 

When there is only a ferroelectric insulator a,n$ np'jo^e, Cd = Cf. 


Chen developed a MFSFET model field effect parameter A. This 

model was only derived for a saturated hysteresis eurve. These models are helpful in 
understanding the theoretical process for the operation of an FFET but do not produce 
data that is consistent with real world FFETs actual measurements. To improve on this 
Bailey and Ho developed a model that partitions the ferroelectric layer into segments. 
The motivation for breaking the ferroelectric layer into segments comes from the fact that 
each section of the layer is under different eonditions. For example, if the gate-source 
voltage is greater the flatband voltage and the drain-source voltage is greater than zero, 
the gate to channel voltage would be constant for the entire length of the channel but the 
semiconductor surface potential would very from Tso at the source to Tsl at the drain. 
This causes the electric field applied to the.ferrpelectric material to vary from 

V GS - ^ (3) at the source to V GS the drain, where d? is the 

ferroelectric thickness. This shows that each segment ofthe ferroelectric material sees a 
different electric field. This indicates that - the pdidrization for each segment of the 
ferroelectric material varies with position. Complicating the model is the fact that the 
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polarization depends not only on the present, epndiiidhs but also on past conditions. This 
approach was presented by Tsividis [4] with respect to ferroelectric polarization. 

The standard Threshold voltage equation musibe modified to account for the charges due 
to polarization of the ferroelectric material.' •>Ifjs|giyeh' by 
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Where Vj is the transistor threshold voltage, V#b is the flatband voltage, Ob -20f + 6 O t 
, where Ob is the voltage drop across the depletion region, Of is the Fermi voltage and O t 
is the thermal voltage, y is the Body effect coefficient, Vsb is the source to body voltage, 
Q’f is the ferroelectric polarization per unit area, and C’f is the ferroelectric capacitance 
per unit area. The difficult part of this equation is determining the overall ferroelectric 
polarization of the ferroelectric material. The model used in this paper divides the length 
of the channel up in to partitions and calculates the polarization numerically given the 
previous polarization, material characteristics, ' arid $ 1 $ voltages present at the transistor 
terminals. 

Figure 1 shows the layout of the model uri$ep conditions. The ferroelectric 

material is divided into segments. For this paperi'tlie number of segments used was 100. 
Increasing the number of partitions did not significantly change the model’s output. 
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Fig. 1. Partitioned MFSFET under operation conditions. 


To obtain the FFET threshold voltage the total amount of charge at the channel interface 
must be determined. This is determined by finding the polarization for each segment 
using ■■ -•<; I';- ■ 
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then summing all of the partitions. The key parameter to solve this equation is the 
surface potential at each partition. The relationship between the channel position and the 
surface potential is given by : v fr -s >v . 

(7) 
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Where x is the channel position, L is the Igri^l?^^ pliannel and /is a function given by 

/(y s )= - v F , + «>,)/ J ( 8 ) 

Since/ (Ts) is nonlinear in Ts, the solution i/dbtaiiied by a numerical root solving 
method of the following equation 
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Once the threshold voltage is determined fo’r the' particular FFET conditions the drain 
current is calculated by 
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Where, . ^ 

V ds is the modified drain to source voltage, .¥b|‘is the gate to source voltage, Vt is the 
transistor threshold voltage, k is an empirical scaling constant, and p, is the electron 
mobility. 


B. Velocity Saturation Model 

The effect of velocity saturation has a very significant affect when the dimension of the 
transistor becomes small. Both the channel length and the oxide thickness have an effect. 
The model uses equations developed by Hodges |2 l The first effect is that of varying the 
effective electron mobility. For high electric fields across the channel oxide a large 
number of carriers are induced in the inversiqiS layef Changing the effective mobility. 

The equation for the effective mobility is given by 
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Where 

© and q are empirical values (3.6 x 10 6 and 1.85 where used). 



The drain current for the linear region is given by ; * 
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Where ■; i: - v VT>V v 

E c is the electric field across the channel and L is the channel length. 


The drain current for the saturation region is given by 



(15) 



= Wv 


sat 



(VgszJtrf 

rLM- ;■ 4 


Where 

V sa t is the saturation velocity. 


Another effect that becomes important when transistor dimensions become very small is 
the direct tunneling as described by Wong When the oxide becomes very t h i n then 
the reduced barrier and increase electric field cause an increase in electrons directly 
tunneling from the gate to the channel. This reduces the desirable characteristics of the 
transistor. The tunneling current density is given by * 
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Where S is the sub-threshold swing, v t is the, thermal voltage, Vt is the threshold voltage, 
tpb and (pbo are 3.1, m ox is 0.4, and a is 0.6, These parameters are for SiC >2 and not for 
PZT. It is assume the PZT parameters would be similar enough to produce similar 
results. 


IV. Results 

A. Verification of model with measured data 

The first result that must be verified is that the new model accurately characterizes an 
actual measured FFET. Measurements from a PZT FFET were compared to the model’s 
predictions. The FFET had a channel length approximately 40 microns long and a 
ferroelectric thickness of about 200 nanometers. The operating gate voltage was +/- 8 
volts. These parameters were input into the model and the resulting output was generated 
in figure 2 for the active mode (Gate voltage GN)iai vo. . 




Fig 2. Model results for drain current for an .FFET in Active mode with L=40,000 

nm and Tf=200 nm. 


This prediction compares quite well to the measured data from the FFET in figure 3. 



Fig. 3. Measured Drain current of actual FFET in Active mode with 1-40,000 nm 

and Tf~200 nm. 


The model also generates data for the remnart mode (Gate voltage OFF after applying 
the gate voltage shown on the X axis). This data is shown in figure 4. 





Fig. 4. Model results for drain current for an FFET in Remnant mode with L=40,000 

nm and Tf=200 nm. 


This again compares favorably to the measured data in figure 5. The remnant measured 
data is affected by the decay of the polarization of the ferroelectric material during the 
time it took to make the measurements. The model does not have any decay time 
modeled between gate pulses. This is why the Upper portion of the measured data seems 
to be lower that the model predicts. 



and Tf~200 nm. 




B. Combined Model Results 

The model was run in both the active and remnant modes for several combinations of 
channel length and ferroelectric film thicknesses. The simulation starts at 40,000 
nanometer channel length because that is the length of the measured FFET. Each 
subsequent simulation has the 1/1 0 th the previous channel length. Channel lengths of less 
than 4 nanometers produced no current in the model. Similarly the ferroelectric film 
thickness varied between 200 nanometers and 3 nanometers. At thicknesses less than 3 
nanometers, the model became unstable and produced unreliable results. Figure 6 shows 
the active mode drain current predicted by the model for a variety of combinations 
channel length and film thicknesses. 



Fig. 6. Simulated drain current for Active mode (Gate voltage = +/- 8V). 
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Similarly, figure 7 shows the remnant mode drain currents predicted by the model for the 
same combinations. 




Fig. 7. Simulated drain current for Remnant mode (Gate voltage = +/- 8V). 


Because the empirical data was taken from a large PZT FFET the model uses gate 
voltages of +/- 8 volts. For future nanoscale transistor different materials and lower 
voltages will be appropriate. Therefore, the simulation was also run with gate voltages of 
+/- 3 volts. Because the PZT does not fully polarize at these voltages, the shapes of the 
drain current versus gate voltage plots have a different shape. These results are shown in 
figures 8 and 9. 



Fig. 8. Simulated drain current for Active mode (Gate voltage = +/- 3V). 




Fig. 9. Simulated drain current for Remnant mode (Gate voltage = +/- 3V). 


C. Direct Tunneling Current 

When the ferroelectric thickness gets very small the lower barrier and the high electric 
fields cause the direct tunneling current to rise to a level that could interfere with the 
normal operation of the transistor. 




Fig. 10. 


Direct tunneling predicted for a lOOnm x lOOnm FFET 


V. Conclusion 

Model simulation suggests that very small ferroelectric transistors are possible. At some 
point the transistors become inoperable due to gate current tunneling into the channel and 
drain current reduction due to reduction in charge of the ferroelectric film, velocity 
saturation and carrier mobility reduction. The point at which these effects become 
overwhelming is about 4 nanometer channel length and 3 nanometer ferroelectric film 
thicknesses. These two parameters are independent but the model becomes unstable if 
either parameter goes below these thresholds. The authors of this paper tried to include 
all major effects on nano-ferroelectric transistors, but because no characterizations of 
actual nano-ferroelectric transistors have be$n made, there may be other effects that 
render the transistor inoperable that are not modeled in this simulation. Future 
simulations need to take more effects into consideration to produce more accurate results. 
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